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Effective theories for the thermal Wilson line are con- 
structed in an SU(N) gauge theory at nonzero temperature. 
I propose that the order of the deconfining phase transition 
for Z(N) Wilson lines is governed by the behavior of SU(N) 
Wilson lines. In a mean field theory, the free energy in the 
deconfined phase is controlled by the condensate for Z(N) 
Wilson lines. Numerical simulations on the lattice, and the 
mean field theory for Z(3) Wilson lines, suggest that about 
any finite temperature transition in QCD, the dominant cor- 
relation length increases by a large, uniform factor, of order 
five. 



A new phase of matter, the Quark-Gluon Plasma, 
might be produced in the collisions of large nuclei at very 
high energies. By asymptotic freedom, the pressure ap- 
proaches the ideal gas value in the limit of high tempera- 
ture, and so it is natural to think of the high temperature 
phase of QCD as a gas of quasiparticles 

It is known, however, that the high temperature phase 
of a purely glue theory is like the low temperature phase 
of a spin system. The magnetization in the high tem- 
perature phase of a SU(N) gauge theory is a Z(N) spin, 
proportional to the trace of the thermal Wilson line |J . 

In this paper I construct effective lagrangians for the 
thermal Wilson line, considered as a full SU(N) matrix, 
as well as its trace. This leads to novel sigma models of 
adjoint SU(N) fields. Although the critical behavior is 
inexorably governed by the fixed point of Z(N) spins j|, 
the SU(N) spins can be important. In particular, they 
help explain why the order of the deconfining transition 
appears to change with N: from second order for N = 2 
|§-§, to weakly first order for N = 3 0-0, to first order 
for N > 4 ]Il| , ^2[ . Further, the picture of the high tem- 
perature phase is turned on its head: the pressure isn't 
due to quasiparticles but is a potential for a con- 

densate of Z(N) Wilson lines. A mean field theory then 
suggests that because the deconfining transition in pure 
glue SU(3) is weakly first order, QCD is near a critical 
point. About the transition, the dominant correlation 
lengths increase by a large factor, of order five ||. 

I concentrate on the pure glue theory; later I argue why 
this is legitimate, using the lattice data and the effective 
theory. The thermal Wilson line is p|,|l3|] 



L(x) = "Pexp ig 



l/T 



Aq(x, t) dr 



(1) 



where V is path ordering, g is the gauge coupling con- 
stant, Aq is the time component of the vector potential 
in the fundamental representation, x is the coordinate for 
three spatial dimensions, and r that for euclidean time 
at a temperature T. The Wilson line in ([!]) is a prod- 
uct of SU (N) matrices, and so is itself a SU (N) matrix, 
satisfying 



L f (x)L(x) = 1 , det(L(x)) = 1 



(2) 



Without quarks, the allowed gauge transformations are 
periodic up to an element of a global Z(N) symmetry [9: 



L(x) 



exp 



(3) 



Q(x) = £l(x, 0). L(x) transforms as an adjoint field under 
local SU (N) gauge transformations in three dimensions, 
and as a vector under global Z(N) transformations. 

Effective theories for L(x) are dictated by the symme- 
tries of (||) . I begin with the nonlinear form, where (||) are 
taken as constraints on L. I construct an effective theory 
in three spatial dimensions, valid for distances 3> l/T, by 
coupling the gauge potentials for static magnetic fields, 
the Ai(x)'s, to the Wilson line, L(x): 



Co 



T 2 aitr|AL| 2 



9* 



(4) 



The first term is the standard lagrangian for static A4 
fields (by choice, A^ has dimensions of mass; all la- 
grangians have dimensions of (mass) 4 ). In the second 
term, I start with electric part of the gauge lagrangian, 
~ tr|Z?iv4o| 2 , and assume that it transmutes into a gauge 
invariant kinetic term for L(x). This is the continuum 
form of the lattice model of Banks and Ukawa jl4j . 

Notice the factor of T 2 in front of the kinetic term 
for L. This arises because the Wilson line is a phase in 
color space, and so every element is a dimensionless pure 
number. Thus in any effective lagrangian, dimensions 
can only be made up by powers of the temperature T. 

Consider the somewhat peculiar limit in which one 
drops the coupling to the Aj's, by taking g — > 0, but 
retains L 7^ 1. Then (^) reduces a nonlinear sigma 
model in three dimensions, with lagrangian ~ tr|<9iL| 2 . 
With the constraints of (^) , the theory is invariant under 
L — > fiiLr^, where fl± and are independent, constant 
SU(N) matrices. This is an enhanced global symmetry 
of SU(N) x SU(N) (times the usual global Z{N) sym- 
metry). As a sigma model, it is possible to impose other 
constraints upon L beyond those of <^). For example, 
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requiring trL to be some fixed number produces a sigma 
model on a symmetric space | fl5f . 

At nonzero coupling, L is simply an adjoint field under 
the local SU (N) symmetry. Even with the constraints of 
(|J) , the reduced symmetry implies that many more terms 
arise: instead of trL being fixed, as for a symmetric space, 
arbitrary traces, such as trL p for integer p, are allowed. 
Mathematically, trL p is related to the trace of the Wilson 
line in higher representations Jl6[ . 

At one loop order, the terms up to fourth order in 
A have been computed. The quadratic term, ~ tr(Ag) 
is the Debye mass for the gluon. For N > 4, 
there are two independent quartic terms, ~ (tr(Ag)) 2 
and ~ tr(Ag) |l9| , p0| ], which represent a potential for Aq. 
From (0), it is easy to turn a potential for Aq into one 
for L: 

& =T 4 ( C2 |trL| 2 + c 4 |trL 2 | 2 + 4(|trL| 2 ) 2 ) . (5) 

Expanding to ~ A\ fixes c 2 = -(4 + 3/tt 2 )/9, c 4 = +(1 + 
3/7r 2 )/36, and c 4 = 0. The rational terms in c 2 and 
c 4 are from the Debye mass, while those ~ l/7r 2 arise 
from quartic terms in the potential for A®. Only two 
constants, c 2 and c 4 , are needed to fit three terms in 
the Aq potential. If Nf flavors of massless quarks are 
included, c 2 and c 4 change, while c 4 is then nonzero. 

The signs of c 2 and c 4 are interesting. As L ~ — g 2 A%, 
a positive Debye mass corresponds to negative c 2 . The 
coupling c 4 is like the quartic term in Aq, and so positive. 
Negative c 2 favors condensation in a direction in which 
| trL | 2 is maximized. This happens when L is an element 
of the center j2lj , 

(L) = exp(2mj/N) £ 1 , (6) 

j = Q...(N — 1). Different j are the usual N degenerate 
vacua of the broken Z(N) global symmetry. 

In (|^) I introduce an expectation value, £q = {£), where 
£ is defined in (|J). In perturbation theory, £q = 1, but 
£q is a function of temperature; it vanishes at the critical 
temperature, T c , and in the confined phase, for T < T c . 

I assume that in the deconfined phase, T > T c , the 
stable vacuum is that which maximizes |trL| 2 , so that L 
condenses as in (^). An expectation value for a field in 
the fundamental representation always breaks the gauge 
symmetry, but uniquely for an adjoint field, a vacuum 
expectation value proportional to the unit matrix does 
not: (^J) is invariant under arbitrary local gauge rota- 
tions. Similarly, the adjoint covariant derivative in (Q) is 
Dih = diL — ig[Ai, so with (^J), the static magnetic 
gluons do not acquire a mass when L condenses, £o ^ 0. 
Thus (^|) is the nonperturbative statement that electric 
screening does not generate screening for static magnetic 
fields gf|(],f§. 

The terms in (|5j) are invariant under a global symmetry 
of U(l). There are also terms which reduce this U(l) to 
Z(N). For N — 3, the simplest examples include 

dot L + c.c. , (trL) 3 + c.c. , trL (trL 2 ) + c.c. . (7) 



The first term, detL, is SU(3) x SU(3) symmetric, while 
the others are only SU (3) symmetric. 

There are also a wide variety of kinetic terms possi- 
ble. These include |9itrL| 2 , <9,;trL 2 | 2 , and |trL| 2 tr|Z?iL| 2 , 
amongst others. At one loop order, the kinetic term in 
(0) is renormalized, and terms such as these may be gen- 
erated; present calculations cannot distinguish [ p3| . Even 
for g — 0, none of these new kinetic terms are invariant 
under SU(N) x SU(N). 

The potential in (ph is only illustrative. In perturba- 
tion theory, one expands about L ~ 1, which does not 
allow one to uniquely fix the coefficients of a potential for 
L. Through numerical simulations, effective theories for 
A [ p0[ , and those for L, could be matched by comparing 
physical correlation lengths at an intermediate tempera- 
ture scale, say at several times the critical temperature. 

Now consider a point of second order transition, where 
£o(T) — > 0. Then powers of L are suppressed, and it is 
sensible to construct a linear sigma model. This is done 
by introducing an "block spin" L, formed by a gauge 
covariant average of L over some region of space [ p4[ . 
Any SU (N) matrix can be written as 

L(x) = £(x)l + 2i£ a (x)t a , (8) 

where t a are the generators of SU(N), a = 1 . . . (A^ 2 — 1). 
For general N, £ and £ a are complex valued, and (||) 
imposes A 2 + 1 constraints. 

I start with the case of two colors, which is special. 
Four constraints of (|^) are satisfied in an especially simple 
manner: the imaginary parts of £ and £ a vanish. This 
leaves one constraint, which is £ 2 + £ 2 a = 1; thus £ and £ a 
form a vector representation of SU(2) x SU(2) = 0(4). 
After averaging, the constraint on the 0(4) norm is lost, 
as is typical in a linear model. Averaging still leaves £ 
and £ a as real valued fields, though. Up to quartic order, 
the most general lagrangian is 

C = \tv{G%) + \{ddf + tr|ZV? (9) 

-mi (£ 2 +£ 2 a )- m 2 £ 2 + Aj (£ 2 + £ 2 a ) 2 + \ 2 £ 4 + A 3 £ 2 £\ . 

The -^-field is a color singlet, while £ a is an adjoint SU(2) 
field. The terms ~ mi and Ai are 0(4) symmetric; with 
the kinetic terms, they correspond to the gauged nonlin- 
ear sigma model of (^|). The other terms, ~ m 2 , A 2 and 
A3, correspond to the potential for the Wilson line in (j|). 
A factor of T has been absorbed into the definition of £ 
and £ a - 

When N > 3, the constraints of (||) are nonlinear. 
Since a sum of two special unitary matrices is not nec- 
essarily special unitary, the average L must be taken to 
be a complex N x N matrix. Thus L includes a complex 
valued, color singlet field, £, which I call a Z(N) spin, 
and a complex valued, color adjoint field, £ a , which I call 
a SU(N) spin. 
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Linear models like @ can be written down for N > 3, 
although there is a plethora of terms. At quartic order 
there is one term which is 0(2iV 2 ) symmetric, (trL^L) 2 , 
another which is SU(N) x SU(N) symmetric, tr(L t L) 2 , 
and terms which are only invariant under SU (N) , such 

as {I = lata) 

{\£\ 2 ) 2 , £tr{ft) 2 £ + c.c. , £ 2 tr(F) 2 + c.c. , 
{txMf , |tr?| 2 , tx{Mf , tr(ff) 2 ? . (10) 

These models give a qualitative picture of the deconfin- 
ing phase transition: I assume that while only the Z(N) 
£-spins condense, (£) = £q ^ 0, that the transition is 
driven by the behavior of the SU (N) £ a -spins. This pic- 
ture is based on the nonlinear model: at weak coupling, 
(|J) dominates other terms, such as (@), by ~ 1/g 2 - Now 
certainly all coupling constants change with T, as can 
be seen from the temperature dependence of the Debye 
mass Nevertheless, I assume that the SU(N) £ a - 

spins dominate right down to the point of the deconfining 
phase transition. The only purpose of terms such as (|5|) 
is to ensure that condensation which respects the local 
SU(N) symmetry, (^), is favored. 

For two colors, the influence of the SU(2) £ a -spins on 
the Z(2) ^-spins is subtle. Assume that only the 0(4) 
symmetric mass, mi, changes. The phase transition in a 
gauged SU (2) model is known from lattice studies of the 
electroweak phase transition P5[ . I assume that one is 
always in an extreme "type-II" regime, so that the second 
order 0(4) transition of the model with g = (the point 
£?2 of fig. (1) in p5[) is washed out by confinement of 
nonabelian gauge fields. The only transition is a point 
at which the Z(2) €-spins become massless; the SU(2) 
£ a -spins are always massive. Lattice studies confirm a 
second order transition in the Z(2) universality class ||. 

For N > 3, the SU(N) £ a -spins can have first order 
transitions. This is because in the absence of gauge fields, 
SU(N) x SU(N) spin models have first order transitions 
both in mean field theory p(i| and in an expansion about 
4 — e dimensions |p5| . As suggested in pj| , in the extreme 
type-II regime, confinement of the gauge fields need not 
wash out the first order transition of SU(N) x SU(N) 
spins (above the point B 3 in fig. (2) of p5[), an d so the 
deconfining transition can remain first order. In particu- 
lar, the transition can be of first order as N — > oo. This 
is in accord with a lattice analysis of Gocksch and Neri 
|P7| , P8| , and contrary to previous speculation [|ll|]3C| l. 

For three colors, this implies that the deconfining tran- 
sition is of first order not only because of cubic invariants 
H, as in (0), but because of the dynamics of SU(3) la- 
spins. Relative to the ideal gas, the latent heat for three 
colors is ~ 1/3 0. As could have been guessed from 
the lattice data alone, perhaps the deconfining transition 
is weakly first order for N — 3 because it is near the 
second order transition for N = 2. Thus it is of value 
to know how the latent heat for N = 4 compares to 
that for N = 3: is it more strongly first order, such as 



~ (N — 2)/N as N — > 2, or more weakly first order, like 
- l/N as N -> oo 0? 

Whatever the order of the deconfining phase transition, 
one can write a mean field theory in which the free energy 
in the deconfined phase is controlled by a potential for 
the Z(N) Wilson lines. For three colors, this is |p|pj}|: 

V = (-2b 2 \£\ 2 + b 3 (£ 3 + (£*f) + (\£\ 2 ) 2 ) b 4 T 4 . (11) 

£ is complex valued, so when b 3 ^ 0, the global symmetry 
is reduced from 0(2) to Z(3). The coupling b 3 must be 
small for the transition to be weakly first order [j3l[ , so 
for now I ignore it, considering the potential just as a 
function of b 2 and 64. This is similar to the case of two 
colors, where £ is a real field, and the potential is just a 
sum of two terms, ~ b 2 £ 2 and ~ 64 £ 4 |32] 1. 

In speaking of the Wilson line, implicitly I assume that 
it is possible to extract a renormalized value |3^] from the 
bare quantity Q. If so, then given £o(T) and the pres- 
sure, one could fit to a potential like (|ll|); for example, 
is it necessary to include higher powers of £ in V? 

The novel aspect of ([ll]) is my insistence that because 
I is a dimensionless field, the dimensions in V must be 
made up by the temperature, T. In mean field theory, 
64 is taken as constant, and 6 2 varies with temperature, 
vanishing at T c . The pressure is given by the minimum 
of the potential, p — b^b^ 4 , and vanishes in the con- 
fined phase, T < T c . That is, with the overall T 4 in the 
potential, the pressure is like a gas of quasiparticles, al- 
beit with a variable number of degrees of freedom, which 
vanish at T c . 

At high temperature, 62 — * 1 so that Iq — ► 1. The 
quartic coupling is fixed by the ideal gas limit: if — 
p/T 4 as T — > 00, 64 = tIqq . Lattice simulations Q [nj 
find that the p/T 4 is relatively flat down to a scale which 
is several times T c , call it nT c ; the same is found from 
resummations of perturbation theory 121. (k might be 
defined as the lowest value where £q w 1 . j Hence I assume 
that 62 and 64 are slowly varying down to kT c . 

Between kT c and T c , I assume that 64 is essentially 
constant, while 62 varies. In particular, the trace of the 
energy momentum tensor, divided by T 4 , is (e— 3p)/T 4 = 
Td(b?,b4)/dT . Lattice simulations find that this quantity 
has a peak just above T c P-|lO||: this is then due to the 
rapid variation of 62 with temperature [j32[ . 

In the Z(N) mean field theory, the pressure includes 
only the contribution of the potential, and nothing from 
fluctuations in the effective fields, either from the Z(N) £- 
spins or the SU(N) ^ a -spins. Fluctuations in these fields 
do, of course, contribute to the pressure at all temper- 
atures. Since by construction the pressure in the mean 
field approximation vanishes for T < T c , one condition 
for its validity is that the pressure in the confined phase 
is small. This is what present lattice simulations find. 
Physically, these fields don't contribute much to the pres- 
sure because they are heavy: the SU(N) £ a -spins always 
so, and the Z(N) ^-spins usually so. For two or three 
colors, the ^-spins do become light in a narrow band in 
temperature about T c , where mean field theory fails [p2| . 
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Ignoring fluctuations about the mean field theory is 
also justified from the viewpoint of an expansion in a 
large number of colors, N — * oo |27],|3(|[35|] . The free 
energy in the confined phase is of order one, while it is 
~ TV 2 in the deconfined phase. The term ~ Y 2 in the free 
energy is due entirely to the condensate, taking b 2 ~ 1 
and &4 ~ Y 2 . Even though there are ~ Y 2 of them, the 
SU(N) ^ a -spins only contribute to the pressure at ~ 1, 
since they are bound into color singlet glueballs. The 
Z(N) ^-spins also contribute ~ 1 to the free energy. 

I have concentrated on the pure glue theory because 
numerical simulations have demonstrated the following 
remar kable property [|§. If p/in^T 4 ) is plotted versus 
T/T c , the resulting curve is nearly universal, and looks 
very similar whether or not there are dynamical quarks 
present. The present model predicts that the pressure is 
the same because the (renormalized fp3| ) Wilson line is 
the same. In terms of the potential, (|ll|), the differences 
in the ideal gas values, n^, are absorbed into 64, with 
the same b 2 {T/T c ). 

Quarks act like a background magnetic field for the 
real part of I [|l4,^6). Because the pure glue transition 
is weakly first order, it is not difficult for quarks to wipe 
out the deconfining transition altogether, leaving either 
a chiral transition, or just crossover behavior. Even so, 
what is relevant here is that for three colors and two or 
three flavors of quarks, the pressure for T < T c is always 
much smaller than that for T > T c ; that is, up, down, 
and strange quarks act like a weak magnetic field for the 
Z{3) ^-spins. 

I thus come to the central physical point of this paper. 
The lattice tells us that the deconfining transition in pure 
glue St/ (3) theory is close to the second order transition 
for SU(2); further, that the effects of quarks are small, 
except close to T c . I suggest that what is important is 
not whether the weakly first order transition persists with 
quarks, but that the nearly second order transition very 
well might. In the pure glue theory, as T — > T c + the ratio 
of the screening mass to the temperature decreases by a 
factor of ten: from - 2.5 at T ~ 2T C , to ~ .25 at T ~ T+ 
Similarly, the string tension at T ~ T~ is ten times 
smaller than that at zero temperature With quarks, 
the increase in the correlation length for £ is presumably 
less, maybe not ten, but perhaps a factor of five or so. 
And most importantly if the pressure below T c is small, 
it might be justified to use the Z{3) mean field theory. 

If the chiral order parameter is then it couples to 
Z{3) £-spins through the coupling +|^| 2 tr . Lat- 

tice simulations find that the chiral and deconfining tran- 
sitions occur at approximately the same temperature. 
This naturally results if this coupling constant is positive, 
as condensation in one field tends to suppress condensa- 
tion in the other. Coherent oscillations in the £-field cou- 
ple to light mesons through such a term, and can produce 
large fluctuations in the average pion momentum ]37| . 

This uniform increase in correlation lengths near T c 
is a unique prediction of the Z(3) mean field theory. 
In quasiparticlc models of the quark-gluon plasma, the 



pressure is tuned to vanish at T c by the introduction of 
a bag constant. In order for the energy to decrease as 
T — ► T c + , though, the quasiparticles must become heav- 
ier, not lighter; that is, instead of increasing, most cor- 
relation lengths decrease Jl],E| . 

At nonzero quark chemical potential, y,, presumably 
there is little change if the quarks are hot and dilute: for 
small jLt/T, the Z{3) £-spins should still exhibit nearly 
second order behavior. I contrast this with the (possible) 
critical endpoint of the chiral transition in the [i—T plane 
p8[ . The correlation length of the sigma meson truly 
diverges at the critical endpoint, but this only occurs 
at one special value of fi. Moreover, the sigma meson 
does not dominate the free energy, nor generic particle 
production. For cold, dense quark matter, fi 3> T, I do 
not see why Z{3) £-spins should dominate the free energy. 

I conclude by noting that the generalization of the De- 
bye mass term, ~ ttvlg, to real scattering processes pro- 
duces hard thermal loops |^] . This is then the first term 
in an infinite series of such terms, continuing ~ trAg, 
etc. The natural expansion is not in powers of Aq, but 
in powers of the Wilson line, as in (||). It is then of great 
interest to know the analytic continuation of the Wilson 
line to real scattering processes EJ . 

I benefited from discussions with K. Eskola, F. Gelis, 
K. Kajantie, F. Karsch, C. P. Korthals-Altes, M. Laine, 
D. Miller, S. Ohta, A. Peshier, D. H. Rischke, M. Wingate 
(for the argument in [0), L. Yaffe, and, especially, M. 
Creutz. This work is supported by DOE grant DE-AC02- 
76CH00016. 
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Define null vectors K — (i, k) and K* — (— i, k), 
and gauge "potentials", 



k 2 = 1, 



A a = 



D ■ K 



1 



D-K* 



With 



(N + N f /2)g 2 T 2 /3 the Debye mass 
squared, and J dflj. the normalized integral over k, 
the usual lagrangian for hard thermal loops is Chti = 
] . This does not look like the expan- 



m 2 J df^tr(A, 
sion of a Wilson line, though. However, an equivalent 
form is 



jChti = m a 



d% tr (.4 A,-K). 



which does look like an expansion of exponentials of A-K, 
and A* ■ K. 
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